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ABSTRACT 

We present deep spectroscopic data for a 24-object subsample of our full 41-object 
z ~ 0.5 radio galaxy sample in order to investigate the evolution of the Fundamental 
Plane of radio galaxies. We find that the low-luminosity, FRI-type, radio galaxies in our 
sample are consistent with the local Fundamental Plane of radio galaxies defined by 
Bettoni et al. when corrected for simple passive evolution of their stellar populations. 
However, we find that the higher luminosity, FRII-type radio galaxies are inconsistent 
with the local Fundamental Plane if only passive evolution is considered, and find 
evidence for a rotation in the Fundamental Plane at z ~ 0.5 when compared with 
the local relation. We show that neither passive evolution, nor a mass-dependent 
evolution in the mass-to-light ratio, nor an evolution in the size of the host galaxies 
can, by themselves, plausibly explain the observed tilt. However, we suggest that some 
combination of all three effects, with size evolution as the dominant factor, may be 
sufficient to explain the difference between the planes. 

We also find evidence for a correlation between host galaxy velocity dispersion 
and radio luminosity at the 97% significance level within our sub-sample, although 
further observations are required in order to determine whether this is different for 
the FRI and FRII radio sources. Assuming that the Mbh — & relation still holds at 
z ~ 0.5, this implies that radio luminosity scales with black hole mass, in agreement 
with previous studies. 

Key words: galaxies: active - galaxies: fundamental parameters - galaxies: nuclei. 



1 INTRODUCTION 

A number of tight correlations exist between the central 
black hole and the host galaxy of an active galactic nucleus 
(AGN): for example, the relation between black-hole mass 
and stellar velocity dispersion (Mbh-ct) (Gebhardt et al. 
2000; Ferrarese & Merritt 2000; Merritt & Ferrarese 2001; 
Giiltekin et al. 2009) and the relation between black hole 
mass and bulge luminosity (e.g. McLure & Dunlop 2001; 
McLure & Dunlop 2002; Giiltekin et al. 2009). The existence 
of these relationships suggests that AGNs have an important 
role to play in the formation and evolution of galaxies via 
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AGN-driven feedback. Various mechanisms have been sug- 
gested for this feedback, including a scenario proposed by 
Silk & Rees (1998) in which quasar outflows limit black hole 
masses, dependent on the depth of the potential wells of 
dark matter halos. 

AGN-driven feedback is now incorporated into models 
of galaxy formation and evolution: for example, the semi- 
analytic model of Croton et al. (2006), where AGN-driven 
feedback is separated into "quasar-mode" (associated with 
the efficient accretion of cold gas) and "radio-mode" (asso- 
ciated with the less efficient accretion of warm gas) feed- 
back. Sijacki et al. (2007) also separate AGN feedback into 
two modes in their full hydrodynamical model: the "quasar 
regime" (corresponding to central black holes with high ac- 
cretion rates) and mechanical feedback (corresponding to 
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central black holes with low accretion rates). Despite the 
nomenclature, radio galaxies can exhibit either "quasar- 
mode" or "radio-mode" /mechanical feedback, dependent on 
the accretion onto the central black hole. There is increas- 
ing evidence (e.g. Hardcastle, Evans & Croston 2007, 2009) 
that the accretion mode of the AGN has a direct link with 
the emission line classification scheme introduced by Hine & 
Longair (1979), where radio galaxies are classified as either 
high-excitation (HEGs) or low-excitation (LEGs) . By study- 
ing the host galaxies of radio-loud AGN spanning luminosi- 
ties from the radio mode (typically lower radio luminosities) 
through to the quasar mode (typically higher radio lumi- 
nosities) we should be able to place important constraints 
on the validity of the feedback mechanisms in the models of 
galaxy formation and evolution. Furthermore, radio galaxies 
reside in the most massive galaxies at all epochs (e.g. Jarvis 
et al. 2001; Willott et al. 2003) and as such allow us to trace 
the evolution of the most massive galaxies across virtually 
all of cosmic time. 

One method of doing this is to study the evolution of 
the Fundamental Plane of radio galaxies. The concept of the 
Fundamental Plane (Djorgovski & Davis 1987; Dressier et 
al. 1987) describes the observation that a three-dimensional 
representation of log R e versus log a versus log I e (where R e 
is the effective radius, a the velocity dispersion and 7 e the 
effective surface brightness of the galaxy) for local ellipti- 
cal and lenticular galaxies yields a plane. The Fundamen- 
tal Plane has a lower scatter than the Faber- Jackson L-a 
relation (Faber & Jackson 1976), and thus was initially of 
interest for its use as an alternative distance determination 
method. However, it was also realised that the fundamen- 
tal plane scalings and small scatter were important for con- 
straining models of elliptical galaxy formation. We may de- 
scribe the fundamental plane as follows: 

R e <x a"/-' 3 . (1) 

Use of the virial theorem with constant mass-to-light 
ratio permits us to reproduce this relation with coefficients 
of a — 2 and j3 — 1 (Faber et al. 1987). However, the ob- 
served fundamental plane does not follow this virial scal- 
ing. For example, j0rgensen, Franx & Kjaergaard (1996) find 
a = 1.24 ± 0.07 and /3 = 0.82 ± 0.02. This "tilt" in the fun- 
damental plane suggests that the mass-to-light ratio (M/L) 
varies as a function of galaxy mass (Robertson et al. 2006) . 
Faber et al. (1987) also note that deviations from the fun- 
damental plane can be caused by other changes in M/L, for 
example due to stellar metallicity or age, or the distribution 
of dark and baryonic matter. Robertson et al. (2006) com- 
ment that each of these effects could, in principle, introduce 
a systematic tilt into the fundamental plane if they vary as 
a function of galaxy mass. J0rgensen et al. (1996) also find 
that the slope of the Fundamental Plane is independent of 
cluster properties such as cluster velocity dispersion and gas 
temperature. 

1.1 The Fundamental Plane of Radio Galaxies & 
AGN 

Bettoni et al. (2001) collected photometric and dynamical 
data for 73 low redshift (z < 0.2) radio galaxies from a vari- 
ety of samples. They find that these local active galaxies lie 



on the same Fundamental Plane as the inactive population, 
with the implication that the gas accretion of the central 
black hole does not influence the global properties of early- 
type galaxies. Wolf & Sheinis (2008) investigate the location 
on the Fundamental Plane of ten nearby (z < 0.3) luminous 
quasars. They find that the radio-loud objects in their sam- 
ple fall on the upper extreme region of the Fundamental 
Plane and have host galaxy properties similar to those of gi- 
ant early-type galaxies, whilst their radio-quiet objects are 
located in the region of the Fundamental Plane occupied by 
normal early-type galaxies and have host galaxy properties 
more similar to intermediate-mass galaxies. 

1.2 Redshift Evolution of the Fundamental Plane 

Work has also been undertaken to examine the redshift evo- 
lution of the Fundamental Plane. The intermediate redshift 
Fundamental Plane (z ~ 0.3-0.4) for cluster (Dokkum & 
Franx 1996; Kelson et al. 2000) and field (Treu et al. 2001) 
early-type galaxies is found to be similar to that observed 
locally. Treu et al. (2001) find an offset (from the local re- 
lation) in their intermediate redshift Fundamental Plane of 
early-type field galaxies. This offset increases with redshift 
and is attributed to an increase in effective brightness of 
these galaxies with respect to a local sample. Woo et al. 
(2004) study the Fundamental Plane of the host galaxies of 
15 AGN out to a redshift of z ~ 0.34. They find that the 
host galaxies of these intermediate redshift BL Lac objects 
and radio galaxies lie on the same Fundamental Plane as 
normal local early type galaxies and the local radio galaxy 
sample of Bettoni et al. (2001), although their higher red- 
shift objects possess an offset in the same fashion as the one 
obeserved by Treu et al. (2001). Woo et al. (2004) also ob- 
serve an increase in the mass-to-light ratio of their galaxies 
of ~40% since z ~ 0.3, also in agreement with the findings 
for normal early- type galaxies. 

At higher redshifts (z ~ 0.8-1.3) various authors (di 
Serego Alighieri et al. 2005; Treu et al. 2005a; J0rgensen et 
al. 2006; J0rgensen et al. 2007; Fritz & Ziegler 2009; Fritz et 
al. 2009; Fritz, J0rgensen & Schiavon 2010) find not only an 
offset but also a rotation in the Fundamental Plane when 
compared to the local relation. For example, J0rgensen et 
al. (2007) find a slope described by a = 0.60 ± 0.22 and 
/3 — 0.70 ± 0.06 for galaxies in their two z ~ 0.8 clus- 
ters whilst di Serego Alighieri et al. (2005) find a slope 
described by a = 0.88 ± 0.16 and j3 = 0.63 ± 0.04 for 
their 0.88 < z < 1.3 early-type galaxies. The latter au- 
thors show at the 90% confidence level that the Funda- 
mental Plane rotates with redshift. This is interpreted as a 
mass-dependent evolution of the mass-to-light ratio, where 
the evolution is faster for less massive galaxies. Fritz et al. 
(2010) study field-type early galaxies and demonstrate that 
the rotation of the Fundamental Plane appears to be inde- 
pendent of the environment of the galaxies. There is thus 
an emerging picture in which variation in the mass-to-light 
ratio causes a redshift-dependent offset in the Fundamental 
Plane (see also van Dokkum & Stanford 2003; van der Wei 
et al. 2005). This variation is mass- dependent, giving rise to 
the redshift-dependent rotation observed in the Fundamen- 
tal Plane. In this picture it is galaxy mass, rather than envi- 
ronment, which plays the major role in determining galaxy 
evolution (Treu et al. 2005b) , and less massive distant galax- 
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ies have much younger stellar populations than more massive 
galaxies [consistent with the downsizing scenario of Cowie 
et al. (1996) in which the most massive galaxies form earli- 
est]. It should, however, be noted that there is not universal 
agreement on this picture: Gebhardt et al. (2003) study the 
Fundamental Plane of 36 field galaxies (21 early-type and 
15 disc galaxies) in the redshift range 0.3-1.0 (with a median 
redshift z = 0.8) and find no difference in the slope of the 
plane when compared to the local relation. 

Robertson et al. (2006) employ a multiphase interstellar 
medium model to investigate the contribution made by var- 
ious different factors to the tilt of the fundamental plane. 
They find that merging gas-rich disc galaxy models (with 
gas fractions f gas > 0.3) including the effects of gas dissipa- 
tion (cooling, star formation and supernova feedback) yields 
a scaling of the Fundamental Plane which is similar to the 
observed infrared Fundamental Plane of Pahre, Djorgovski 
& de Carvalho (1998). Changing the redshift of the progen- 
itor galaxies within the range z = 0-6 has little effect on the 
plane produced. The inclusion of gas dissipation has the ef- 
fect of decreasing the total-to-stellar mass ratio (M tota ;/M*) 
in the central regions of the galaxies. This decrease varies as 
a function of galaxy mass (with a greater decrease in lower 
mass systems) thus producing the required tilt in the Fun- 
damental Plane. Re-merging the remnants of these gas-rich 
galaxy mergers leaves the Fundamental Plane largely intact. 
Including the effects of feedback from accreting supermas- 
sive black holes in their models makes little difference to 
the Fundamental Plane inhabited by the merger remnants. 
However, use of these full-model simulations allow Robert- 
son et al. (2005) to reproduce the power-law scaling of the 
Mbh-o relation between z — 6 and z — 0. The remnant 
remergers of Robertson et al. (2006) preserve the Mbh-o 
relation (albeit with increased scatter). 



1.3 The Fundamental Plane at z ~ 0.5 

In this paper we present spectroscopic data for a sample of 
radio galaxies at z ~ 0.5. We use a direct spectral fitting 
procedure to obtain the velocity dispersions of our galaxies. 
The effective radii and surface brightnesses for our sample 
have already been determined and published (see McLure et 
al. 2004). Combining our velocity dispersion results with this 
data thus permits us to study the Fundamental Plane at z ~ 
0.5 and draw comparisons with the well documented local 
relation. Our results also allow us to investigate any possible 
correlation between the low-frequency radio luminosity at 
151 MHz, L151, and the velocity dispersion. This enables 
us to investigate (via the Mbh-o relation) the possible link 
between radio luminosity and black hole mass in radio-loud 
AGN (e.g. Nelson & Whittle 1996; McLure & Jarvis 2004; 
Jarvis & McLure 2006). 

We begin in Section [2] of the paper by describing our 
sample, before detailing our observations and data reduction 
in Section [3] Section [4] describes our spectral fitting and 
aperture correction techniques, followed by a description of 
our Fundamental Plane fitting in Section [5] We present and 
discuss our findings in Section|6]before concluding in Section 
[7] Unless otherwise stated, throughout the paper we assume 
a standard cosmology in which Ho = 70 km s -1 , Qm = 0.3 
and Qa = 0.7. 
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Figure 1. The logarithm of the 151-MHz radio luminosity, L151, 
versus redshift for the ZP5 (z ~ 0.5) objects in this paper (filled 
symbols). We also use open symbols to show the local radio galax- 
ies from Bettoni et al. (2001; 2 ~ 0.03). In each case the sym- 
bols correspond to the 151-MHz radio luminosities. Objects with 
L151 < 10 25 ' 3 W Hz -1 sr — 1 are shown as stars. Objects with L151 
between 10 25 ' 3 and 10 25 ' 8 W Hz -1 sr — 1 are shown as squares 
whilst objects with L151 between 10 25,8 and 10 26 3 W Hz -1 sr" 1 
are shown as inverted triangles. Objects with L151 between 10 26 ' 3 
and 10 26 ' 8 W Hz" 1 sr~ 1 are shown as triangles, and objects with 
^151 > 10 26 ' 8 W Hz" 1 sr" 1 are shown as circles. We do not 
show five of the objects of Bettoni et al. (0131-360, 0325+024, 
PKS 2322-122, 0112+084 and ginll6) for which we have been 
unable to locate published radio luminosity data. 



2 THE SAMPLE 

Our full 41-object z ~ 0.5 radio-galaxy sample (the ZP5 
sample) consists of all of the narrow-line radio galaxies in 
the redshift interval 0.4 < z < 0.6 from four complete, 
low-frequency selected radio surveys; 3CRR (Laing, Riley 
& Longair 1983), 6CE (Eales et al. 1997; Rawlings, Eales & 
Lacy 2001), 7CRS (Lacy et al. 1999; Willott et al. 2003) and 
TexOx-1000 ("TOOT"; Hill & Rawlings 2003; Vardoulaki et 
al. 2010). The choice of 0.4 < z < 0.6 was motivated by sev- 
eral factors. In particular, a sample at z ~ 0.5 provides a 
full three decades in radio luminosity (including some of the 
most radio luminous objects in existence). It therefore allows 
us to construct subsamples of objects, divided by radio lumi- 
nosity, containing sufficient numbers of objects to investigate 
correlations of various properties with radio luminosity. In 
addition, the use of 0.4 < z < 0.6 ensures that satisfactory 
data can be obtained (both Hubble Space Telescope imaging 
and spectroscopic data) without excessively long integration 
times. Full details of the sample can be found in McLure et 
al. (2004) and Herbert et al. (2010). 

In this paper we use a 24-object subsample of the full 
ZP5 sample (see Section [3] and Table [T]). Spectroscopic ob- 
servations were made using the William Herschel Telescope 
(WHT) or the Gemini North telescope. Objects were chosen 
at random from the full sample, aside from the constraint of 
2 < 0.5 for objects observed with the WHT. In Figure [T] we 
show the radio luminosity-redshift plane for the ZP5 objects 
in this paper. 
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Table 1. The z ~ 0.5 sample. Column 1 lists the radio galaxy names and columns 2 and 3 list the J2000 source coordinates. Column 4 
lists the object redshifts, column 5 the logarithm of the 151-MHz luminosities in units of WHz _1 sr _1 and column 6 the effective radii in 
kiloparsccs from McLure et al. (2004). The effective radii are derived from fits in which the Sersic parameter is allowed to vary. Column 
7 lists the mean ij-band surface brightness within the effective radii in mag arcsec 2 , also from McLure et al. (2004). Column 8 lists the 
aperture corrected velocity dispersion of our best fit and column 9 the reduced \ 2 of this fit, whilst column 10 gives the telescope used 
for the spectroscopic observations. 
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3 OBSERVATIONS AND DATA REDUCTION 
3.1 WHT Spectroscopy 

In August 2002 we used the ISIS spectrograph on the 4.2- 
metre William Herschel Telescope to obtain spectra for 15 
of the objects in the full sample. The objects were chosen to 
have z < 0.5 so that satisfactory signal to noise ratios could 
be obtained without excessively long integration times. This 
constraint aside, the objects were chosen at random from the 
objects in the full sample that were visible on the dates of 
observation. Details of the objects observed with the WHT 
can be found in Table [1] We used the R316R grating on the 
red arm together with a 1.0 arcsec slit resulting in a resolu- 
tion of 3.7 A FWHM. This corresponds to a velocity reso- 
lution of FWHM = 185 km s" 1 at 6000A which results in a 
minimum measurable velocity dispersion of a ~ 80 km s . 
The central wavelength used was dependent on the redshift 
such that the spectra contain the Mgfe complex to deter- 
mine the velocity dispersion. Each object was observed for 
a total integration time of 7200s, split into three 2400s ob- 
servations to allow the removal of cosmic rays. The individ- 



ual spectra were bias subtracted, illumination and fiat-field 
corrected, and wavelength and flux-calibrated using stan- 
dard IRAF routines. Coaddition was done using an average 
sigma-clipping routine and the galaxy was extracted using 
an aperture defined at the FWHM of the galaxy light profile 
and at full-width zero intensity to ensure both high signal- 
to-noise ratio in the former and total light spectra for the 
latter. 



3.2 GMOS Spectroscopy 

Between August 2008 and June 2009 we used the Gemini 
Multi-Object Spectrograph (GMOS) on the 8.1-metre Gem- 
ini North telescope to obtain spectra for nine of the objects 
in the full sample. Objects with z > 0.5 were permitted 
and were chosen at random (given RA constraints) from the 
objects in the ZP5 sample left unobserved after the WHT 
observations. Details of the objects observed using Gemini 
can be found in Tabled] The B600_G5303 grating (and sub- 
sequently the new B600_G5307 grating) was used together 
with a 0.75 arcsec slit resulting in a spectral resolution of 
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4.0 A FWHM. This corresponds to a velocity resolution of 
FWHM = 190 km s" 1 at 6200A which results in a min- 
imum measureable velocity dispersion of a ~ 80 km s _1 . 
Within each observation block, two 600s observations were 
made with a central wavelength of 610nm, and a further two 
600s observations were made using a central wavelength of 
605nm. The use of two different central wavelengths allows 
us to cover the gaps between the CCDs. Each observation 
block was repeated four times per source. Thus each object 
was observed for 9600s in total. The CMOS spectra were 
reduced (bias subtraction, cosmic ray subtraction, flat field 
correction, wavelength calibration, sky subtraction, flux cal- 
ibration and spectrum extraction) using the gmos tools in 
the gemini IRAF package. 



An advantage of explicity looping over the velocity dis- 
persion is that it allows us to obtain the standard errors on 
a from the \ 2 distribution. We convert the \ 2 matrix fx?) 
(corresponding to the variation of \ 2 with a for the best- 
fitting template for each galaxy) into a probability distribu- 
tion matrix (Pi) and normalise. We then inspect the proba- 
bility distribution function in order to determine the range 
in velocity dispersion which contains 68.27% of the probabil- 
ity distribution function and thus obtain the standard errors 
on a. As a check we note that this procedure results in the 
same standard errors as produced by A\ 2 = 1. We note that 
three objects (3C341, 6C1303+3756 and TOOT1630+4534) 
have reduced x 2 values > 1.5. The errors on these objects 
are therefore undoubtedly underestimated (see Table [T}. 



4 SPECTRAL FITTING 
4.1 Our Fitting Procedure 

In order to determine the velocity dispersion of each galaxy 
we perform a \ 2 minimisation, fitting each galaxy spectrum 
to a set of templates. The template set comprises 23 Elodie 
stellar templates Q and (following Greene & Ho 2006) 34 
G and K stars from the old open cluster M67 observed 
with SDSS. The templates are broadened using a series of 
Gaussians whose width depends on the velocity dispersion 
(a) and wavelength. We exclude from our fit regions of the 
galaxy spectra that feature emission lines ionised by the cen- 
tral engine, as well as skylines. The effect of broadening due 
to the instrumental resolution is accounted for in our fitting 
procedure, and the effect is typically within the uncertainties 
of our fitted velocity dispersions. 

To implement our fitting we fit each galaxy in turn, us- 
ing all 57 templates and with the velocity dispersion incre- 
menting in steps of 1 km s . At each step we minimise the 
reduced x 2 °f the fit using the Amoeba routine in IDL. We 
then inspect the resulting matrix of reduced \ 2 values in or- 
der to find the minimum reduced \ 2 an d thus the best fitting 
template and velocity dispersion. Using simulated spectra 
generated from our stellar template library, we find that our 
minimum measureable velocity dispersion is < 100 km s , 
well below the expected velocity dispersions of our radio 
galaxy hosts. We are unable to determine reliable measure- 
ments of the velocity dispersion for the cases where the mean 
signal-to-noise ratio per resolution element (measured be- 
tween rest frame 4400A and 4600A where spectral features, 
in particular the Mg6 complex, dominate the fitting) is less 
than ~10. On this basis we exclude our fits to 3C16, 3C427.1, 
6C0825+3407, 6C0850+3747 and 7C0219+3423. Since this 
is based merely on the quality of the data we do not expect 
this to introduce any notable selection effects. In addition, 
we are unable to obtain a sensible fit in the case of 3C457 
due to a combination of the signal-to-noise ratio and emis- 
sion lines from the AGN. Given that this is the only object 
lacking a fit due to the AGN, we are confident that this 
has a negligible effect on our results. Our final fits to the 
remaining objects are shown in Figure [21 Spectra for which 
we were unable to determine a reliable measurement of the 
velocity dispersion are shown in the Appendix. 

1 http://atlas.obs-hp.fr/elodie/ 



4.2 Aperture Correction 

The galaxies in our z ~ 0.5 sample are smaller on the sky 
than those of a local sample. Consequently the aperture size 
for our sample is of order 20 times greater than the aperture 
size for a local sample. This has the effect of including more 
of the galaxy in the slit, resulting in a lower determined ve- 
locity dispersion than would be the case for a local sample. 
Therefore, in order to be able to compare our results with 
existing results for local samples we perform an aperture cor- 
rection on our derived velocity dispersion values and their 
errors. We perform the aperture correction using the pro- 
cedure described by j0rgensen, Franx & Kjaergaard (1995), 
where our values of a are corrected to a circular aperture 
with an aperture diameter of 1.19A kpc (equivalent to 3.4 
arcsec projected onto a galaxy in the Coma cluster). The 
average correction is 7.1%. 

Our derived velocity dispersion values, including aper- 
ture correction, can be found in Table [T] Table [T] also con- 
tains the McLure et al. (2004) values for the major-axis ef- 
fective radii, r e , and the mean .R-band surface brightnesses 
within r e , (fJ,)e, of the galaxies in the sample. The values of 
(/i) e have been corrected for passive evolution as described 
in McLure et al. (2004). 



5 FUNDAMENTAL PLANE FITTING 
5.1 Initial Plane Fitting 

In order to fit a plane to our data we adopt the formulation 
of the fundamental plane used by Robertson et al. (2006), 
namely 

R e oc a a I~ fi , (2) 

where R e is the effective radius and a the velocity dis- 
persion as before. I e is the mean surface brightness within 
R e , in units of L /pc 2 , and is related to our previous quan- 
tity (/j,} e : 

log/ e = -0.4(( M ) e -fc), (3) 

where for the R-beaid k — 26.4 (see j0rgensen et al. 
1996) . In order to fit the fundamental plane we rewrite Equa- 
tion [2] in the following form: 
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Figure 2. Our final fits. The left hand panels show the galaxy spectrum (top) and the fitted template (bottom). A downward shift of 0.3 
has been applied to the template spectra. The shaded areas indicate regions of the spectra removed from the fitting due to the presence 
of emission (rest frame 3859-3880A, 3963-3978A, 4095-4112A, 4332-4378A, 4674-4698A and 4854-4875A) or telluric absorption (5570- 
5585A, 6864-6911A, 7166-7323A and 7590-7694A) features. The additional masked regions in 3C172, 6C1303+3756 and 7C1731+6638 
were due to artefacts on the chip. The right hand panels show the \ 2 minimisation as a function of a Aperture corrections have been 
applied as described in the text. 
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Figure [2] continued. 
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Figure [2] continued. 
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log R e = a log a - /3 log J e - 7, (4) 

where 7 is a constant. For a plane described by Equation 
[4] using coordinates (log r e , log a, log I e ) the residual 1Z of 
point (pi,P2,J>3) perpendicular to the plane is given by: 

U = ~ (-pi + ap 2 - Pp 3 - 7) , (5) 

where 

S = yjl + a 2 + /3 2 . (6) 

We note that 12 of our 18 fitted objects have sym- 
metric errors on logo - . To fit a plane we therefore ini- 
tially adopt symmetric errors in logr e , logo - , and logic. 
For the six sources which do not have a symmetric er- 
ror on log a we use the mean of (log (a + a upP er) — log a) 
and (logcr — log (a — ai ower )), where a upper and (flower COr 
respond to the upper and lower errors on a listed in Table [T] 
The values of r e and I e taken from McLure et al. (2004) and 
used here have linear errors of 10%. We use the same proce- 
dure as above to estimate symmetric errors in log space for 
these quantities. In order to fit a plane to our results we use 
a technique similar to the one used by di Serego Alighieri et 



al. (2005), where we use a finely sampled grid over a, /3 and 
7. At each point we use Equations [5] and [6] to calculate the 
residual perpendicular to the plane of each of our objects. 
We use our symmetric errors in log space in order to de- 
termine the error on each residual. We define a \ 2 statistic 
thus: 



x 2 = E^. m 

i a Ki 

where the summation is over the objects we are fitting 
the plane to, IZi is the residual of object i perpendicular to 
the plane, and an i is the standard error on the residual. 
Finally we find the values of a, j3 and 7 which minimise this 
statistic and thus determine the best fitting plane. Our best 
fit fundamental plane can be found in Table [2] 

j0rgensen et al. (1996) and Bettoni et al. (2001) fit a 
plane by minimising the sum of the absolute residuals. How- 
ever, by using our x 2 -based technique, we can determine the 
standard errors on a, /3 and 7 in a similar fashion to that 
employed in Section f4.1l We convert our three-dimensional 
X 2 matrix into a probability distribution matrix and nor- 
malise. We then determine the volume that contains 68.27% 
of the probability distribution function and from this the 
standard errors on a, /3 and 7. The reduced-x 2 also provides 
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Table 2. Values of a, and 7 for a fundamental plane fitting 
the ZP5 data and described by the equation logr e = a log a — 
ft log I e — 7. The first row gives the results in the case where we 
assume symmetric errors in log space (Section 15,11 1. The second 
row gives the results from the full Monte Carlo (MC) fitting us- 
ing squared residuals (Section 15,211 . The same procedure, except 
weighting the minimised residuals by the error, produces the re- 
sults listed in the third row. The results obtained by minimising 
the sum of the square residuals and the sum of the absolute resid- 
uals without the Monte Carlo technique are given in the fourth 
and fifth rows respectively (Section l5,2l l. 
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us with a measure of the goodness of our fit. Alongside these 
advantages we do, however, emphasize that this technique 
does not make use of the exact errors on r e , a and I e and 
that small approximations are required. However, we show 
in Section 15.21 that the plane we obtain using this method 
is consistent with the plane produced using a method which 
involves no approximations. 



5.2 Monte Carlo Simulations 

In order to fit a plane (with errors) to our data, using the ac- 
tual error distributions on a, we employ a Monte Carlo tech- 
nique. For each of our objects we randomly sample 12,000 
values from each of the r e , I e and a probability distribu- 
tion functions (PDFs). We use sufficient samples to produce 
smooth final distributions in a, j3 and 7. The procedure out- 
lined in section [4~T1 provides the a PDF. For each sample set 
we again grid over a range of a, /3 and 7 and obtain the best 
fitting values of these parameters for 12,000 simulations by 
minimising the sum of the squared residuals according to 
equation [5] We obtain the standard errors by finding the 
range containing 68.27% of the PDF for each parameter. 
Table [2] lists our fitting results. 

Next we repeat this procedure but rather than minimise 
just the sum of the square residuals, we minimise the sum of 
the square residuals weighted by the square of the error on 
the residual (where we estimate the residual error using the 
procedure described in Section l5.1|) . Finally, we fit a plane to 
the data by minimising the sum of the squared residuals, and 
another by minimising the sum of the absolute residuals (as 
per j0rgensen et al. 1996; Bettoni et al. 2001), but without 
any Monte Carlo sampling in either case. The results from 
all of these plane fittings can be found in Tabled and in all 
cases our Fundamental Plane parameters are consistent. 



5.3 Local Plane Fitting 

In order to compare our 2 ~ 0.5 plane with results for local 
active galaxies we use the data presented in Bettoni et al. 
(2001). This comprises both new data and data from the 
literature on their sample of low redshift radio galaxies. The 



median redshift of their sample is 0.03. We note that we use 
the same procedure as Bettoni et al. to apply the aperture 
corrections and thus our results can be compared directly. 
The first row of Table [3] contains the published plane fitting 
of Bettoni et al. (2001), converted from from their plane pa- 
rameterisation to ours using Equation [3] The second row of 
Table[3jgives our fit to the data of Bettoni et al. (2001) using 
their fitting procedure (minimising the sum of the absolute 
residuals perpendicular to the plane). We exclude from our 
data set the objects 0053-016 and 0431-134 from Ledlow & 
Owen (1995), for which we have been unable to locate pub- 
lished redshifts. 

In order to compare our plane fitting with the data of 
Bettoni et al. (2001) we correct the scalelengths in the data 
from their cosmology (Ho = 50 km s _1 , Qm = 1 and Q.a 
— 0) to ours. Once again we exclude the objects 0053-016 
and 0431-134 from Ledlow & Owen (1995) due to the lack of 
redshift information. We fit the Fundamental Plane to the 
corrected data by minimising the sum of the absolute resid- 
uals and then repeat, this time minimising the sum of the 
squared residuals (see Table [3]). Finally, we use the proce- 
dure described in Section [5.11 to fit the Fundamental Plane 
to the corrected data of Bettoni et al. (2001) using a % 2 fit- 
ting statistic, yielding uncertainties on the plane fitting as 
well as an indication of the goodness of fit. Uncertainties on 
a have been published for the data presented by Bettoni et 
al. (2001) and Smith, Heckman & Illingworth (1990). For 
the remainder of the data we use a typical linear error of 
10%. We also assume linear errors of 10% on r e and I e . The 
fitting results can be found in Table [3] 



6 DISCUSSION 

6.1 The Local Fundamental Plane 

In Figure [3] we show a two-dimensional projection of the 
Fundamental Plane that we fit to the cosmologically cor- 
rected data of Bettoni et al. (2001) using their fitting pro- 
cedure (row three of Table [3]). As well as showing the data 
of Bettoni et al. in this projection we also show our ZP5 
objects with symbols corresponding to their radio luminosi- 
ties. The y-axis uncertainties for the ZP5 objects were deter- 
mined by assuming errors as in Section 15.11 and combining 
in quadrature according to the equation of the plane. It can 
be seen in Figure [3] that our objects with radio luminosity 
L151 < 10 25 ' 3 W Hz -1 sr" 1 appear to lie on the local Funda- 
mental Plane. The traditional division between lower radio 
luminosity Fanaroff-Riley type I (FRI) and higher radio lu- 
minosity Fanaroff-Riley type II (FRII) radio galaxies falls 
at L151 = 10 25 ' 3 W Hz" 1 sr" 1 (Fanaroff & Riley 1974), and 
this is also where there appears to be a difference in the evo- 
lution (Clewley & Jarvis 2004; Sadler et al. 2007), and thus 
these objects are probably FRI type radio galaxies, although 
some overlap in the populations is inevitable (see Herbert et 
al. 2010). We note that the majority of the sample of Bet- 
toni et al. (~90% of the objects which have published radio 
luminosities) also have I/151 < 10 25 ' 3 W Hz -1 sr _1 and are 
thus classified as FRI type radio galaxies. We therefore find 
evidence that the FRI type objects in our sample can in- 
habit the local Fundamental Plane of similar objects simply 
by passively evolving from z ~ 0.5 to the present day. 
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Table 3. Values of a, j3 and 7 for a fundamental plane fitting the data of Bettoni et al. (2001) and described by the equation logr e = 
a logo - — /SlogTe — 7. The first row lists the plane fitting results published by Bettoni et al. (2001), where we convert from their plane 
parameterisation to ours using Equation [3] The second row lists our fit to their data obtained by minimising the sum of the absolute 
residuals (excluding the two objects discussed in the text), and the fitting described in the third row is the same but using the data of 
Bettoni et al. (2001) corrected into our cosmology. Row four gives the results of our fit to the corrected data by minimising the sum of 
the squared residuals, and the fifth row gives the results where we assume symmetric errors in log space and minimise the \ 2 statistic. 
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Figure 3. A 2D projection of the Fundamental Plane as fitted to the cosmologically corrected data of Bettoni et al. (z ~ 0.03) by 
minimising the sum of the absolute residuals perpendicular to the plane. The data of Bettoni et al. is shown by open symbols. We use 
filled symbols to show the ZP5 data (z ~ 0.5), corrected for passive evolution and aperture corrected as described in the text. In each case 
the symbols correspond to the 151-MHz radio luminosities. Objects with L151 < 10 25 ' 3 W Hz -1 sr — 1 are shown as stars. Objects with 
£151 between 10 25 3 and 10 25 8 W Hz -1 sr — 1 are shown as squares whilst objects with L151 between 10 25 8 and 10 26 ' 3 W Hz -1 sr _1 
are shown as inverted triangles. Objects with L151 between 10 26 ' 3 and 10 26 ' 8 W Hz -1 sr _1 are shown as triangles, and objects with 
^151 > 10 26 s W Hz -1 sr _1 are shown as circles. Objects from Bettoni et al. for which we have been unable to locate a published radio 
luminosity are shown as crosses. We show the typical assumed error on the data of Bettoni et al. in the top left-hand corner (see Section 

E3). 
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Many of the remaining (FRII type, L151 > 10 25 3 
W Hz -1 sr" 1 ) objects also lie close to the plane and within 
the scatter displayed by the data of Bettoni et al. (2001). 
However, we note that our FRII type radio galaxies exhibit 
a tendency to lie above the local Fundamental Plane, some 
by a reasonable distance (although there are also objects 
lying a reasonable distance below the plane). Passive evo- 
lution is not sufficient to explain how these objects would 
move from their position at z ~ 0.5 onto the 2 = Fun- 
damental Plane. The FRII type objects of Bettoni et al. 
are consistent with the plane occupied by the rest of the 
z ~ sample, although there are insufficient objects with 
high radio luminosity to study this properly. Therefore, in 
contrast to our FRI type objects, we find evidence that the 
host galaxies of our z ~ 0.5 FRII type objects must undergo 
evolution above and beyond their passive evolution to move 
them onto the local Fundamental Plane (Section 16. 3[) . We 
find no clear correlation within the FRII type group between 
radio luminosity and the offset from the local Fundamental 
Plane, suggesting that this is indeed an evolutionary effect 
(rather than one related to the radio luminosity). 

We also study the offsets of our ZP5 objects from the 
local Fundamental Plane of Figure [3] as a function of radio 
morphology and spectral classification. Details of our classi- 
fication schemes can be found in Herbert et al. (2010). We 
find that those objects with large offsets are predominantly 
'Classical Double' (FRII) HEG sources, consistent with the 
above result based on radio luminosity and with a similar 
strength of trend. It is also worth emphasizing that we find 
no link between disturbed morphology from a recent merger 
(see McLure et al. 2004) and the size of the offset. 

6.2 The ZP5 Fundamental Plane 

In Figure [4] we show a two-dimensional projection of the 
Fundamental Plane fitting the ZP5 data assuming symmet- 
ric errors in log space (row one of Table [5]) . It is apparent 
that our ZP5 objects do inhabit a z ~ 0.5 Fundamental 
Plane, which is common to both the FRI type and the FRII 
type objects in our sample. 

The reduced-x 2 of the ZP5 plane fitting shown in Figure 
f3] is 2.9 (tabulated, along with other reduced-x 2 values, in 
Table 2| . The reduced-x 2 obtained when we compare the 
ZP5 data to the local Fundamental Plane of Figure [3] (a = 
1.76, ft = 0.713, 7 = 1.5) is 10.2 (whereas comparing this 
local plane to the data of Bettoni et al. 2001 yields a reduced- 
X 2 of 3.4). For comparison we also fit a plane to our ZP5 data 
that follows the Kormendy (fj, e - r e ) relation (Kormendy 
1977) by setting a = and fitting for ft and 7. Our best fit 
in this case (ft — 0.73, 7 = —2.7) has a reduced-^ 2 of 5.1. 
We find that our ZP5 fit is a noticeably better fit to the ZP5 
data than the local Fundamental Plane is. It is interesting 
to note from row one of Table [2] and row three of Table [3] 
that this ZP5 fit is significantly different from the fit to the 
local data. The most prominent difference between the local 
Fundamental Plane and the z ~ 0.5 Fundamental Plane is 
that the a exponent is much smaller for the z ~ 0.5 plane. 
That is to say, the power law index that relates the effective 
radius and the surface brightness to the velocity dispersion 
is far smaller at z ~ 0.5 than in the local universe. The 
velocity dispersion is related to the mass (M) of the galaxy, 
and thus we infer that at z ~ 0.5 the mass and size of these 



Table 4. Reduced-x 2 values. The first row contains the reduced- 
X 2 values yielded by a comparison of the local Fundamental Plane 
(a = 1.76, P = 0.713, 7 = 1.5) to our data and that of Bettoni 
et al. (2001), whilst the second row shows the reduced-x 2 values 
from comparing the ZP5 plane (row one of Table [2]l to the two 
data sets. The third row contains the reduced-x 2 values obtained 
by comparing the size-evolved ZP5 plane to the two data sets 
(Section 16.31 1. where the ZP5 data is also ajusted to account for 
size-evolution. 
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powerful radio galaxies are less closely related than in the 
local universe. However, despite the smaller a exponent in 
our ZP5 relation, a comparison of the the reduced-x 2 values 
for our normal ZP5 fit and our Kormendy fit suggests that 
our z ~ 0.5 Fundamental Plane is indeed a plane involving 
a rather than merely a fi e - r e relation; a fit with a = 
is noticeably worse than one where a dependence on a is 
allowed. 

Inspecting Figure [4] we find that, whilst the objects of 
Bettoni et al. are offset with regard to the ZP5 Fundamental 
Plane, there nevertheless appears to be a good agreement 
between the two data sets. We therefore perform a plane 
fitting (assuming symmetric errors in log space) to the two 
data sets combined. However, the resulting plane (a — 1.23, 
ft — 0.61, 7 = 0.5) yields reduced-x 2 values of 4.8, 6.1 and 
3.5 when compared to the combined data, the ZP5 data and 
the data of Bettoni et al. respectively. We therefore find that 
a plane fitted to the combined data is not acceptable. 

In Section 14.11 we noted that the errors on a for some 
of our objects are undoubtedly underestimated. This con- 
tributes in part to the large reduced-x 2 values of the fits. 
In order to test whether these underestimated errors are 
leading us to reject acceptable plane fittings we repeat our 
fittings and comparisons, but requiring a minimum error of 
15% on r e , cr, and I e . The reduced-x 2 values we obtain from 
comparing the data to the various plane fittings are shown 
in Table(5] We find, with a reduced-x 2 °f 1-3, that the local 
plane (Figure [3} is an acceptable fit to the data of Bettoni 
et al.. We also find that the ZP5 plane (Figure |3| is an ac- 
ceptable fit to the ZP5 data (reduced-x 2 of 1.2). However, 
we find that the local plane is not an acceptable fit to the 
ZP5 data (with a reduced-x 2 of 3.5), and nor is the ZP5 
plane an acceptable fit to the local data (with a reduced-x 2 
of 2.4). Finally, a plane fitted to the two data sets combined 
is also still rejected (reduced-x 2 values of 1.9, 2.6 and 1.3 
when compared to the combined data, the ZP5 data and 
the data of Bettoni et al. respectively) . 

Taken altogether, our results therefore suggest that ra- 
dio galaxies with lower and higher radio luminosities may 
share the same Fundamental Plane, but that there is sub- 
stantial evolution in this plane between z — and z ~ 0.5. 

6.3 Evolution of the Host Galaxies? 

One explanation for the difference between the local Funda- 
mental Plane of Bettoni et al. (2001) and our z ~ 0.5 Fun- 
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Figure 4. A 2D projection of the Fundamental Plane as fitted to the ZP5 data (z ~ 0.5), assuming symmetric errors in log space. We 
show the ZP5 objects and the local radio galaxies from Bettoni et al. (2001; z ~ 0.03) in this plane projection, with symbols as in Figure 
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Table 5. Reduced-x 2 values where the data used in the plane 
fittings and for the comparisons have errors > 15%. The rows 
and columns are as in Table 2] 
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damental Plane is evolution of the host galaxies. We have 
already shown in Sections 16. II and 16. 21 that passive evolution 
is not sufficient to explain the difference between the planes, 
or at least for the high radio-luminosity population. Other 
possibilities include an evolution in the mass-to-light ratio 
or a size evolution of the host galaxies. We discuss each of 
these in turn. 

Various authors have found evidence for a mass- 
dependent evolution of the mass-to- light ratio (Section ll,2l) . 
and this in principle could explain the observed rotation be- 
tween our ZP5 Fundamental Plane and the local one. We 
use our two plane fittings (row one of Tableland row three 
of Table [3| to calculate the mass-dependent evolution of 
the mass-to-light ratio that would be required to transform 



our z ~ 0.5 plane into the local one (assuming r c and a 
remain constant and expressing a in terms of r e and M). 
However, we find an implausibly large dependence on the 
z ~ 0.5 luminosity would be required (greater than 10 or- 
ders of magnitude, with a correspondingly small dependence 
on the mass), and we therefore discount the possibility that 
this, by itself, could explain the difference between the two 
planes. 

In Figure [5] we show histograms of the mass-to-light 
ratio distributions for the ZP5 and Bettoni samples. A two- 
dimensional K-S test on the two distributions rejects at the 
96% level the null hypothesis that the samples are drawn 
from the same distribution. We therefore find tentative evi- 
dence of an evolution in the mass-to-light ratio which, whilst 
not explaining the observed rotation between the z — and 
the 2 ~ 0.5 planes by itself, may make some contribution 
to the tilt. Considering only the FRII type objects in our 
sample we reject the null hypothesis that this sample and 
that of Bettoni et al. are drawn from the same distribution 
at the 98% level. However, we find no clear evidence that the 
distribution of mass-to-light ratios is different for our FRI 
type objects and the sample of Bettoni et al. (rejecting the 
null hypothesis at only the 86% level), consistent with these 
objects lying on the local Fundamental Plane of Figure [3] 
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Figure 5. Histograms of the mass-to-light ratio distributions for 
the ZP5 (z ~ 0.5) and Bettoni (z ~ 0.03) samples. 



We next turn to consider the possibility of size evolution 
of the host galaxies. There exists evidence for the growth 
of early-type galaxies from higher redshifts to the present 
day: for example, the observational work of van der Wei et 
al. (2008) and the model of Hopkins et al. (2010), but see 
Mancini et al. (2010). We use the size evolution model of van 
der Wei et al. (2009), where mergers cause the mass and ef- 
fective radius to increase by the same factor. Thus a remains 
constant (since a oc -y/ M/R), along with the mass-to-light 
ratio, whilst r e and 7 e vary. We characterise the size evolu- 
tion using a factor Sr e where r e (z = 0.5) = Sr e ■ r e (z — 0). 
We investigate what value of Sr e gives the smallest offset 
between the ZP5 plane and the local data by varying Sr e 
(and adjusting r e and I e accordingly), refitting a plane us- 
ing the procedure described in Section [5.11 and minimising 
the reduced-^ 2 of the fit of the local data of Bettoni et al. 
(2001) to this plane. For each value of 5r e we adjust the 
amount of passive evolution such that our ZP5 objects still 
satisfy the Kormendy relation of the bottom panel of figure 
5 in McLure et al. (2004). We note that using the same Sr e 
factor for all galaxies simply produces a uniform translation 
of the objects with respect to the logr e and log/ e axes and 
thus affects the offset but not the slope of the corresponding 
ZP5 plane. However, for each value of Sr e we re-fit a plane 
(rather than simply adjust the 7 plane parameter) in order 
to account for the uncertainties in the a and /3 plane pa- 
rameters. The best fit is obtained with Sr e — 0.39 (although 
the reduced-x 2 is fairly constant between 0.3 and 0.6) and 
is shown in Figure [6] This is in contrast to a value of ~0.91 
expected from the model of van der Wei et al. (where we 
consider only the effects of mergers and exclude the effects 
of late-type galaxies turning into early-types) but the upper 
end of our range is in good agreement with the model of 
Khochfar & Silk (2006) who predict a value of Sr e ~ 0.63- 
0.64 for the most massive galaxies. The size evolution im- 
plied by 5r e = 0.39 would result in a mean scalelength for 
our objects of 37.7 kpc. This is similar to the size of Bright- 
est Cluster Galaxies (BCGs) in the local universe (Graham 
et al. 1996; McLure et al. 2004), which suggests that the host 
galaxies of powerful radio galaxies may evolve into BCGs. 

However, as discussed above, whilst applying a uniform 
size evolution factor to our ZP5 objects can reduce the off- 
set to the local galaxies, it does not affect the rotation. This 
is seen in Figure [6] where the tilt of the plane is still not 



consistent with the local data of Bettoni et al. (2001). In- 
deed, whilst the reduced-x 2 obtained by comparing the size- 
evolved ZP5 data to the plane of Figure E] is 2.9 (Table g]), 
or 1.2 when we use errors > 15%, when we compare the 
data of Bettoni et al. to this plane the reduced-^ 2 is 4.8 
(or 1.9 with > 15% errors). Similarly, fitting a plane to the 
combined data of Bettoni et al. and the size-evolved ZP5 
data yields a plane (a — 1.32, /3 = 0.61, 7 = 0.7) whose 
reduced-x 2 values are 4.5, 5.3 and 3.4 when compared to 
the combined data, the ZP5 data and the data of Bettoni 
et al. respectively. Using errors > 15% these values reduce 
to 1.8, 2.3 and 1.3 respectively. We therefore find that size 
evolution of the host galaxies by a uniform factor is not suf- 
ficient, in and of itself, to explain the difference between the 
z — and the z ~ 0.5 Fundamental Planes, since it doesn't 
address the difference in rotation between these two planes. 

In order for the tilt to be explained by size evolution we 
would require a mass (and size) dependent evolution of the 
host galaxies. Once more comparing the ZP5 plane fitting 
in row one of Table [2] with the local fitting in row three of 
Table and taking 5r e (as defined above) to be a function 
of mass, we find Sr e oc rl' 52 M~ 211 (neglecting any other 
forms of evolution). This seems more plausible than the 
mass-dependent mass-to-light ratio evolution, although the 
actual evolutionary factors required for our objects (0.09- 
86.7, calculated using the same method) are not. A mass- 
dependent size evolution could be driven by environmental 
effects (where the most massive galaxies are found only in 
virialized clusters where mergers are suppressed), which we 
will investigate (Herbert et al., in preparation) using deep 
multi-band imaging data. 

We therefore find that, whilst a mass-dependent size 
evolution may be the dominant effect behind the evolution 
of the Fundamental Plane, it is not sufficient in and of itself 
to explain the observed rotation between the z ~ 0.5 and 
z — planes. We suggest that size evolution may combine 
with passive evolution and a mass-dependent evolution of 
the mass-to-light ratio to produce the observed rotation. 



6.4 A Link with Recent Star Formation? 

We next consider the possibility that our powerful ZP5 radio 
galaxies already inhabit the local Fundamental Plane, but 
that the effects of recent star formation make them appear 
to lie off it. In other words, is it possible that our objects 
that lie below the local Fundamental Plane could have a 
bright extended disc of recent star formation that would 
increase the observed effective radius? Conversely, is it pos- 
sible that our objects that lie above the local plane could 
have a bright nucleocentric region of recent star formation 
that would decrease the observed effective radius? If the star 
formation terminates by z = then the observed effective 
radius would be altered in both cases, possibly moving the 
objects back on to the local relation. 

We test this hypothesis using the strength of the 4000A 
break for our objects; the D„ (4000) indices (taken from Her- 
bert et al. 2010) can be used as a measure of recent star 
formation, where younger populations have smaller 4000A 
breaks, as discussed in Herbert et al. (2010). In Figure [7| we 
show the residual perpendicular to the plane of Figure[3]ver- 
sus D„(4000) for the ZP5 objects. We see that those objects 
that fall significantly below the local Fundamental Plane, as 
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Figure 6. A 2D projection of the Fundamental Plane as fitted to the ZP5 data (z ~ 0.5), size-evolved by a factor Sr e = 0.39. We show 
the size-evolved ZP5 objects and the local radio galaxies from Bettoni et al. (2001; z ~ 0.03) in this plane projection, with symbols as 
in Figure [3] 



well as the two furthest above it, all have lower D n (4000) 
indices and thus younger stellar populations. This lends a 
degree of support to the recent star formation hypothesis, 
although we note that some of the objects that fall furthest 
from the plane have no evidence of recent star formation 
whilst some falling on or near the plane do have evidence 
for younger stellar populations. 

In order to test the hypothesis further we extract the 
spectra for our objects using a range of apertures across 
the galaxy. We calculate D n (4000) indices in each case and 
look for any trends between D n (4000) and extraction radius 
for each object. However, we find no such trends within the 
uncertainties on D n (4000) and the amount of possible con- 
tamination discussed in Herbert et al. (2010). We cannot 
therefore, with our current data, prove star formation in ex- 
tended discs and nucleocentric regions as an explanation for 
why some of our objects lie away from the local Fundamental 
Plane. 

6.5 Correlations with Radio Luminosity 

Evidence of a correlation between black hole mass and radio 
luminosity was first presented by Franceschini, Vercellone 
& Fabian (1998), who showed that both the nuclear and 
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Figure 7. The residual perpendicular to the plane of Figure [3] 
versus D n (4000) for the ZP5 objects (z ~ 0.5). The symbols are 
as used in Figure [3] 

the total radio fluxes of their sample of nearby radio galax- 
ies displayed a remarkably tight dependence on black hole 
mass. This result has been confirmed by more recent work 
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(e.g. McLure et al. 1999; Lacy et al. 2001; McLure & Jarvis 
2004) whilst evidence has also emerged of a connection be- 
tween black hole mass and radio loudness (e.g. Laor 2000; 
McLure & Dunlop 2001; McLure & Jarvis 2004), although 
other studies have found no evidence of such connections 
(e.g. Ho 2002; Woo & Urry 2002). 

In order to investigate the possible link between radio 
luminosity and black hole mass for the objects in our sample 
we calculate the black hole mass from a using the Mbh- 
a relation for elliptical galaxies published by Giiltekin et 
al. (2009) and propogate errors in the standard fashion by 
assuming symmetric errors on log a as above. It is important 
to note that, whilst the majority of the previous studies 
have used the high-frequency (5-GHz) radio luminosity, we 
use the extended low-frequency (151-MHz) radio luminosity. 
This may be significant, since the 151-MHz radio luminosity 
is closely connected with the time-averaged kinetic energy of 
the jets (e.g. Rawlings & Saunders 1991) and the effects of 
beaming are smaller on the 151-MHz radio luminosity than 
on the 5-GHz radio luminosity (e.g. Jarvis & McLure 2002). 

In Figure [5] we show L151 versus black hole mass. Fig- 
ure [8] provides evidence in support of a link between radio 
luminosity and black hole mass: using the Spearman rank 
correlation coefficient we find a positive correlation between 
the L151 radio luminosity and Mbh at a significance level of 
97%. This is in good agreement with the L\*,\-Mbh relation 
for the whole ZP5 sample found at a significance level of 97% 
by McLure et al. (2004) where the black hole masses were es- 
timated via the MsH-MBuige relation of McLure & Dunlop 
(2002). However, in contrast to McLure et al. (2004), who 
find that the significance of the correlation is increased to 
99.5% by excluding the TOOT objects, we find that exclud- 
ing the TOOT objects decreases the correlation significance 
to 84%. However, this could be purely due to the decrease in 
the number of objects we use for our analysis; 24 compared 
to the full sample of 41 used in the work of McLure et al. 
(2004). Additional spectroscopy on the remaining 27 radio 
galaxies would be needed to confirm this. 

We also find a positive correlation between L151 and r e 
(Figure [9]). Applying the Spearman rank correlation coef- 
ficient yields a significance of 98%. It is interesting to note 
that the TOOT objects (triangles) — FRI type objects with 
L151 < 10 25 - 3 W Hz -1 sr" 1 — do not appear to follow this 
relation; for the TOOT objects r e appears to be largely inde- 
pendent of I/151. This is consistent with our findings for the 
Fundamental Plane, where the low-luminosity radio sources 
are consistent with being able to passively evolve onto the 
local relation. 



7 CONCLUSIONS 

In this paper we have presented deep spectroscopic observa- 
tions of 24 powerful radio galaxies at z ~ 0.5. We have used 
a direct fitting procedure to determine the velocity disper- 
sions of our objects and combined these with the host galaxy 
data from McLure et al. (2004) to study the Fundamental 
Plane. 

• We have shown that the FRI type objects in our sample, 
once corrected for passive evolution, lie on the Fundamental 
Plane inhabited by the local radio galaxies of Bettoni et al. 
(2001) . The FRII type objects in our sample do not lie on the 
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local plane. However, whilst the majority of the sample of 
Bettoni et al. are lower radio luminosity (FRI type) objects, 
the FRII type objects in their sample are consistent with 
the local plane. Furthermore, our FRI type objects and our 
FRII type objects appear able to share the same plane at z ~ 
0.5. We therefore suggest that radio galaxies with lower and 
higher radio luminosities may share the same Fundamental 
Plane, but that there is substantial evolution in this plane 
between z — and z ~ 0.5. 

• One explanation for the tilt observed between the z — 
and the z ~ 0.5 Fundamental Planes is evolution of the host 
galaxies. We consider passive evolution, a mass-dependent 
evolution in the mass-to-light ratio and an evolution in the 
size of the host galaxy, but find that none of these effects can, 
by themselves, plausibly explain the observed rotation. We 
suggest, however, that some combination of all three effects, 
with size evolution as the dominant factor, may be sufficient 
to explain the difference between the planes. Indeed, indirect 
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evidence for size evolution for the high-luminosity sources 
comes from the significant correlation (98%) between the 
effective radius and the radio luminosity for these sources, 
whereas the FRI sources do not appear to align with such a 
relation. 

• We also consider the possibility that our radio galax- 
ies do, in fact, lie on the local Fundamental Plane (once 
passive evolution has been corrected for), but that a bright 
extended disc or nucleocentric region of recent star forma- 
tion alters the observed effective radius, thus causing some 
of our objects to appear to lie off the local relation. Using 
the D„(4000) index as an indicator of recent star formation, 
we see some hints that this may be having an effect, but the 
evidence is far from conclusive. 

• We find evidence at the 97% level of a correlation be- 
tween galaxy velocity dispersion and the radio luminosity, 
suggesting that radio luminosity scales with black-hole mass. 
This is in line with previous work on this sample using 
the host galaxy luminosity as a proxy for black-hole mass 
(McLure et al. 2004). Unlike McLure et al., however, we do 
not find that this correlation becomes stronger in the ab- 
sence of the low-luminosity (FRI-type) radio sources. We 
attribute this to the lack of objects in the present study 
(24) compared with the analysis of the full sample of 4f in 
McLure et al.. 

In future work (Herbert et al. in preparation) we will 
use deep multi-band imaging data to investigate the envi- 
ronments of the full 2 ~ 0.5 sample and how they relate to 
the power of the AGNs and the masses and star formation 
histories of the host galaxies. 
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APPENDIX 

In Figure[T0]we show the spectra of the six objects for which 
we were unable to obtain a reliable measurement of the ve- 
locity dispersion (Section l4.ll) . 
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Figure 10. Spectra of the six objects for which we were unable to obtain a reliable measurement of the velocity dispersion. 3C457 was 
rejected due to a combination of the signal- to- noise ratio and emission lines from the AGN. The remainder were rejected on the basis of 
their signal-to- noise ratios alone (see Section 14.11 1. 



